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Abstract 

In  this  paper,  we  address  the  problem  of  consistency  checking  for  Euclidean  spatial  constraints. 

A  dimension  graph  representation  is  proposed  to  maintain  the  Euclidean  spatial  constraints  among 
objects.  The  basic  idea  is  to  project  the  spatial  constraints  on  both  X  and  Y  dimensions,  and  the 
dimension  graph  is  constructed  on  each  dimension.  By  using  the  dimension  graph  representation, 
the  problem  of  consistency  checking  is  then  transformed  to  a  graph  cycle  detection  problem.  The 
consistency  checking  can  be  achieved  with  0(N+E)  time  as  well  as  space  complexity,  where  N  is  the 
number  of  spatial  objects,  and  E  is  the  number  of  spatial  predicates  in  the  constraint.  The  proposed 
approach  to  consistency  checking  for  spatial  constraints  is  faster  than  0(N2)  when  the  number  of 
predicates  is  much  smaller  than  N2  and  there  are  few  disjunctions  in  the  spatial  constraint.  The 
dimension  graph  and  consistency  checking  algorithm  can  be  used  for  points,  intervals  and  polygons 
in  2  dimensional  space.  The  algorithm  can  guarantee  the  global  consistency. 

Keywords:  Euclidean  spatial  constraint,  consistency  checking,  dimension  graph,  directional  relation¬ 
ship 
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1  Introduction 


A  spatial  database  [12,  14,  23]  management  system  aims  at  the  effective  and  efficient  management  of 
data  related  to  a  space  such  as  the  physical  world  (geography,  urban  planning,  astronomy);  parts  of 
living  organisms  (anatomy  of  the  human  body);  engineering  design  (very  large  scale  integrated  circuits, 
the  design  of  an  automobile  or  the  molecular  structure  of  a  pharmaceutical  drug);  and  conceptual 
information  space  (a  multi-dimensional  decision  support  system,  fluid  flow,  or  an  electro-magnetic  field). 
The  distinguishing  features  of  a  spatial  database  management  system  are  the  use  of  complex  data 
types  like  points,  lines  and  polygons  to  represent  spatial  objects  and  the  existence  of  many  potential 
relationships  between  spatial  objects. 

Consistency  checking  is  an  important  concept  to  maintain  the  integrity  of  databases  in  general 
and  spatial  databases  in  particular.  Consistency  checking  is  the  process  of  identifying  contradictory 
information  in  a  database.  For  example,  if  A,  B  and  C  are  three  spatial  objects  and  if  B  is  west  of  A 
and  C  is  west  of  B  and  if  the  database  indicates  that  C  is  east  of  A  then  the  information  is  inconsistent  . 
The  existence  of  many  potential  spatial  relationships  implies  that  consistency  checking  in  the  context  of 
spatial  databases  is  more  challenging  vis-a-vis  its  traditional  relational  counterpart. 

Consistency  checking  can  also  be  used  for  spatial  query  processing  via  semantic  query  optimization[21] 
and  reasoning  for  spatial  qualitative  relationships.  For  example,  given  a  spatial  query  S  and  a  set  of  of 
spatial  constraints  SC,  if  the  query  S  is  inconsistent  with  respect  to  SC  then  the  answer  to  S  is  null. 
Currently  most  of  consistency  checking  is  based  on  Allens ’s  algorithm  [3]  which  was  originally  devised 
for  checking  temporal  relationships  on  one-dimensional  objects.  Spatial  relationships  are  typically  for¬ 
mulated  among  multi-dimensional  objects  and  straightforward  extension  of  Allen  algorithm  for  spatial 
relationships  are  prohibitively  expensive. 

In  this  paper,  we  address  the  problem  of  consistency  checking  for  directional  spatial  constraints  in 
the  two-dimensional  Euclidean  space.  We  propose  a  dimension  graph  representation  for  maintaining  the 
spatial  constraints  among  objects.  Basically,  the  spatial  constraints  is  projected  on  each  dimension(X 
and  Y),  the  constraints  that  must  be  satisfied  for  each  dimension(X/Y)  are  maintained  in  different 
graph(X/Y  graph).  One  graph  records  the  constraints  on  one  dimension.  The  problem  of  constraint 
consistency  checking  is  transformed  to  a  graph  cycle  detection  problem  on  the  dimension  graph.  The 
cycle  detection  could  be  solved  by  traversing  the  graph  in  linear  time.  As  we  will  see  in  later  sections, 
the  proposed  consistency  checking  algorithm  is  efficient  in  terms  of  both  time  and  space.  The  algorithm 
also  guarantees  the  global  consistency  of  system. 

1.1  Spatial  Data  Model 

Recent  reports[12,  14,  23,  1]  have  described  the  accomplishments  of  spatial  database  research  and  have 
prioritized  research  needs.  A  broad  survey  of  spatial  database  requirements  and  an  overview  of  research 
results  is  provided  by  [23,  12,  1,  20]. 

A  spatial  data  model  is  defined  by  geometric  entities,  spatial  operations  on  these  entities  and  spatial 
relationships  between  them.  The  three  basic  geometric  entities  are  point,  line  and  polygon  which  rep- 
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resent  spatial  objects  in  zero,  one  and  two  dimension  respectively.  To  facilitate  rapid  query  processing 
polygons  are  often  represented  by  their  Minimum  Bounding  Rectangles(MBRs).  An  MBR  of  a  spatial 
object  is  the  smallest  axis  parallel  rectangle  which  covers  the  polygon. 

The  spatial  relationships  between  the  geometric  entities  are  categorized  according  to  the  mathemat¬ 
ical  properties  of  the  relationships.  For  example,  topological  relationships  like  overlap  and  adjacent  are 
relationships  which  are  invariant  under  a  change  of  coordinate  or  projection  system.  Other  types  of 
relationships  include  Metric,  Set  and  Directional.  Some  important  relationships  categories  and  their 
representative  example  are  shown  in  1.  In  this  paper  we  will  focus  on  consistency  checking  of  absolute 
directional  relationships  like  North,  South,  East,  West,  Northeast,  Northwest,  Southeast  and  Southwest. 
These  relationships  can  be  defined  between  any  combination  of  point,  line  and  polygon  objects  but  we 
will  restrict  our  attention  to  homogeneous  point,  line  and  MBR  pairs.  For  constraints  among  intervals, 
we  deal  with  Allen’s  13  relationships  [3].  The  constraints  are  represented  in  terms  of  disjunctions  and/or 
conjunctions  of  spatial  predicates.  The  predicates  are  logically  atomic. 


Data  model 

Operator  Group 

Operation 

Vector  Object 

Set-Oriented 

equals,  is  a  member  of,  is  empty,  is  a  subset  of,  is  disjoint 
from,  intersection,  union,  difference,  cardinality 

Topological 

boundary,  interior,  closure,  meets,  overlaps,  is  inside,  covers, 
connected,  components,  extremes,  is  within 

Metric 

distance,  bearing/angle,  length,  area,  perimeter. 

Direction 

east,  north,  left,  above,  between. 

Network 

successors,  ancestors,  connected,  shortest-path 

Table  1:  A  Sample  of  Spatial  Operations 


1.2  Problem  Definition 

In  this  paper,  we  intend  to  explore  the  consistency  checking  for  Euclidean  qualitative  spatial  constraints 
among  points,  intervals,  and  2-D  Minimum  Bounding  Rectangles(MBRs)  objects  in  spatial  database. 

Consistency  Checking  Problem: 

Given:  A  collection  of  0-th  order1  spatial  constraints  in  terms  of  disjunctions  and/or  conjunctions 
of  spatial  predicates 

Find:  Consistency,  i.e. ,  return  TRUE  if  the  constraints  are  consistent,  False  otherwise 
Objective:  Reduce  computation  complexity 
Constraint:  (a)  2-D  Euclidean  space. 

(b)  Two-dimensional  extended  objects  are  approximated  by  MBRs 

(c)  Spatial  objects  refereed  to  given  spatial  constraints  are  of  homogeneous  types, 

(i.e.,  point  pairs,  interval  pairs,  or  MBRs  in  2D) 


10-th  order  spatial  constraints  means  that  there  are  no  free  variables  but  only  constant  objects 
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Consider  an  example  of  direction  constraints  among  point  objects  A,  B,  and  C.  Assume  one  constraint 
says  A  is  Southeast  of  B,  B  is  Northwest  of  C,  and  A  is  Northeast  of  C.  This  constraint  is  consistent, 
and  hence  the  result  will  be  TRUE.  An  example  of  the  possible  spatial  configurations  satisfying  this 
constraint  is  shown  in  Figure  1. 

North  B. 


East  •  ^ 

C 

Figure  1:  One  possible  spatial  configuration  for  point  objects  A,  B,  and  C 

Suppose  we  have  another  set  of  constraint  among  A,  B,  and  C:  A  is  strictly  north  of  B,  B  is  Northwest 
of  C,  and  A  is  Northeast  of  C.  This  constraint  is  inconsistent  since  there  does  not  exist  any  spatial 
configuration  of  A,  B,  and  C  satisfying  this  constraint.  The  consistency  checking  algorithm  should 
return  FALSE. 

1.3  Related  work  and  our  contributions 

Most  of  previous  study  on  consistency  checking  is  based  on  Allen’s  consistency  checking  algorithm  [3]  for 
constraints  among  intervals.  The  basic  approach  for  consistency  checking  is  to  use  transitive  closure  al¬ 
gorithm,  which  incurs  high  order  of  time  and  space  complexity.  Hernandez  [13]  presented  mechanisms  to 
maintain  the  consistency  of  a  knowledge  base  of  spatial  information  based  on  a  qualitative  representation 
of  2-D  positions.  His  approach  improved  Allen’s  algorithm  by  using  the  heuristic  of  rich  structure  of  the 
spatial  domain.  Bowman  and  etc. [6,  5]  addressed  the  problem  of  consistency  checking  between  multiple 
viewpoints  using  strategies  based  on  unification.  Some  other  work[15,  10,  17]  focused  on  the  problem 
of  consistency  checking  for  more  general  constraints,  such  as,  1-th  order  constraints.  Manandhar[16] 
discussed  deterministic  consistency  checking  for  LP  constraints.  Beneventano  and  etc. [4]  focused  on 
the  problem  of  providing  a  theoretical  framework  for  consistency  checking  of  integrity  constraints  in  a 
complex  object  database  environment. 

The  consistency  checking  algorithm  is  dependent  on  the  set  of  spatial  predicates.  This  paper,  we  deal 
with  basic  direction  predicates  for  absolute  directions[22].  There  are  some  other  directional  predicates. 
The  research  work  on  direction  modeling  has  been  carried  out  in  several  areas  such  as  geographic 
information  systems  and  image  analysis.  Most  of  the  studies  is  on  how  to  capture  the  semantics  of 
direction  relations,  and  further,  how  to  do  spatial  reasoning  based  on  direction  predicates  [8,  9,  11].  There 
are  two  major  direction  reference  frames  used  to  model  direction  in  2D  space:  the  cone-based  model[18], 
and  the  projection-based  model[9,  11].  Frank[2]  compared  these  two  models  and  found  the  projection- 
based  reference  frame  to  be  better  in  many  aspects.  The  most  common  way  to  model  directions  between 
extended  objects  is  through  the  object’s  Minimum  Bounding  Rectangle(MBR),  where  direction  relations 
are  obtained  by  applying  Allen’s  [3]  interval  relations  along  the  x  and  y  axis,  in  which  case,  169  different 
relations[8]  can  be  distinguished.  We’ll  briefly  describe  Allen’s  Algorithm  since  it  is  basic  to  many 
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algorithms. 


1.3.1  Allen’s  Propagation  Algorithm 

Allen[3]  summarized  thirteen  mutually  exclusive  relationships  to  express  any  possible  relationships  be¬ 
tween  intervals  as  shown  in  table  2. 


Relationships 

symbol 

Symbol  for  Inverse 

Pictorial  Example 

X  before  Y 

< 

> 

XXX  YYY 

X  equal  Y 

— 

XXX 

YYY 

X  meets  Y 

m 

mi 

XXX YYY 

X  overlaps  Y 

0 

oi 

XXX 

YYY 

X  during  Y 

d 

di 

XXX 

YYYYYY 

X  starts  Y 

s 

si 

XXX 

YYYYY 

X  finishes  Y 

f 

A 

XXX 

YYYYY 

Table  2:  Thirteen  possible  relationships  proposed  by  Allenfl] 


In  Allen’s  work[3],  the  relationships  between  intervals  are  maintained  in  a  network  where  each  node 
Ni  represents  individual  interval  i,  and  each  arc  N(i,j)  is  associated  with  possible  relationships  between 
the  corresponding  interval  pair  i  and  j.  The  basic  algorithm  he  used  for  maintaining  relationships  was 
propagating  new  relationships  by  computing  the  transitive  closure  of  the  relationships  between  intervals. 

Figure  2  shows  the  network  used  by  Allen  for  consistency  checking.  2(a)  is  the  network  for  two 
inputs,  i.e.  S  overlaps  or  meets  L ,  and  S  is  before,  meets,  is  metby,  or  after  R.  After  the  second  input 


Figure  2:  Examples  of  Allen’s  algorithm 

was  added,  the  algorithm  computed  the  constraint  between  L  and  R,  and  the  resulting  network  is  shown 
as  2(b).  If  we  add  a  new  fact  L  overlaps,  starts,  or  is  during  R,  we  need  to  propagate  its  effect  through 
the  network,  and  thus  obtaining  the  resulting  network  2(c). 

As  explained  in  Allen’s  paper[3],  the  time  complexity  of  this  algorithm  is  calculated  as:  13  x 
(\  for  ^  intervals,  i.e.,  0(N'2).  The  space  requirement  for  the  algorithm  is  also  0(N2). 


4 


As  Allen  stated  in  [3],  one  problem  with  this  algorithm  is  that  it  does  not  detect  all  inconsistencies 
in  its  input.  Quote  the  phrases  in  [3]  ”In  fact,  it  only  guarantees  consistency  between  three  node 
subnetworks.  There  are  networks  that  can  be  added  which  appear  consistent  by  viewing  any  three  nodes, 
but  for  which  there  is  no  consistent  overall  labeling  of  the  network.”  In  other  words,  the  algorithm  can 
not  guarantee  global  consistency. 

1.3.2  Our  Approach 

In  this  paper,  we  propose  a  new  strategy  to  process  consistency  checking  for  Euclidean  spatial  constraints 
among  objects.  We  use  a  geometric  approach  by  incorporating  the  spatial  domain  information  in  consis¬ 
tency  checking.  We  propose  dimension  graphs  to  maintain  the  spatial  constraints  among  objects.  Each 
conjunctive  constraint  is  projected  on  both  dimensions,  and  a  dimension  graph  is  constructed  for  the 
constraint  on  each  dimension.  For  spatial  constraints  in  general  format,  the  constraints  can  be  converted 
to  its  Disjunctive  Normal  Form(DNF)[19]  ,  and  dimension  graphs  are  constructed  for  each  conjunction. 
By  using  the  dimension  graph  representation,  the  problem  of  constraint  consistency  checking  is  then 
transformed  to  a  graph  cycle  detection  problem  on  each  dimension  graph.  The  cycle  detection  could  be 
solved  efficiently  with  0(N+E)  time  as  well  as  space  complexity,  where  N  is  the  number  of  spatial  objects, 
and  E  is  the  number  of  spatial  predicates  in  the  constraint.  Recall  that  Allen’s  algorithm  has  a  time 
and  space  complexity  of  0{N2).  The  proposed  approach  to  consistency  checking  for  spatial  constraints 
is  faster  when  the  number  of  predicates  is  much  smaller  than  N 2  and  there  are  few  disjunctions  in  the 
spatial  constraint.  The  dimension  graph  and  consistency  checking  algorithm  can  be  used  for  points,  and 
MBRs  in  2  dimensional  space.  Since  the  algorithm  returns  TRUE  if  and  only  if  the  dimension  graph 
of  at  least  one  conjunction  contains  no  cycle,  which  means  there  exist  at  least  one  consistent  overall 
constraint.  The  algorithm  can  guarantee  the  global  consistency. 

1.4  Scope  and  Outline 

In  this  paper,  we  address  the  problem  of  consistency  checking  for  spatial  constraints  in  the  two- 
dimensional  Euclidean  space.  We  deal  with  only  0-th  order  constraints.  We  focus  on  the  qualitative 
constraints  among  objects,  which  include  topological  and  direction  relationships(See  Table  1).  Distance- 
based  constraints  are  not  discussed.  We  consider  the  constraints  among  point  objects,  interval  objects, 
and  region  objects  approximated  by  MBRs.  We  focus  on  addressing  the  consistency  checking  for  homoge¬ 
neous  types  of  objects,  i.e.,  the  constraints  among  point  objects,  or  intervals,  or  MBRs.  The  consistency 
checking  for  the  constraints  specified  among  mixed  types  of  objects(e.g.  the  constraints  between  a  point 
and  an  interval)  is  out  of  the  scope  of  the  paper,  which  may  be  addressed  in  the  future  work.  We  only 
discuss  a  specific  set  of  predicates  defined  in  the  paper.  The  predicates  are  defined  in  Euclidean  space. 
Different  predicate  sets  or  in  different  space  may  be  different.  Some  constraints  that  are  inconsistent  in 
Euclidean  space  may  be  consistent  in  spherical  space.  The  consistency  checking  for  predicates  in  other 
space  is  out  of  the  scope  of  this  paper. 

The  organization  of  this  paper  is  as  follows:  In  section  2,  we  propose  dimension  graph  representation 
for  the  conjunctive  constraints  among  points,  intervals.  The  consistency  checking  for  conjunctive  con- 
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straint  based  on  the  dimension  graph  representation  is  introduced  in  section  3.  In  section  4,  we  discussed 
the  dimension  graph  construction  and  consistency  checking  for  conjunctive  constraints  among  MBRs. 
Finally,  the  consistency  checking  for  constraints  in  general  format  is  described  in  section  5.  The  paper 
ends  with  conclusions  and  recommendations  for  future  work. 


2  Dimension  Graphs  for  Conjunctive  Spatial  Constraints 

In  this  section,  we  will  describe  the  construction  of  dimension  graph  for  conjunctive  spatial  constraints 
among  points  and  intervals.  The  basic  idea  is  to  projected  the  conjunctive  spatial  constraint  onto  each 
dimension  and  record  the  spatial  constraints  that  much  be  satisfied  on  each  dimension  in  a  dimension 
graph  respectively.  The  dimension  graph  for  points  in  2D  space  contains  X-graph  and  Y-graph.  The 
dimension  graph  for  intervals  in  ID  space  contains  only  one  graph.  We  also  analyze  the  computation 
complexity  for  the  dimension  graph  construction  algorithm. 


2.1  Dimension  Graphs  for  Constraints  Among  Point  Objects 

We  start  by  defining  a  set  of  absolute  direction  predicates  for  point  objects  in  terms  of  coordinates.  Here, 
We  assume  the  global  coordinate  systems  aligned  with  the  reference  frame  of  absolute  directions,  i.e. , 
North  aligns  with  y-axis,  and  East  align  with  x-axis.  The  definition  for  each  predicate  is  given  in  Table 
3.  The  first  column  of  the  table  enumerates  the  direction  predicates.  The  second  and  third  columns 


Direction  predicates 

Ax  5  Bx 

Ay  ,  By 

SP(A,  B) 

= 

= 

North(A,  B) 

= 

> 

South{A,  B) 

= 

< 

Easi{A,  B) 

> 

= 

West(A,B) 

< 

= 

NE(A,  B) 

> 

> 

NW(A,B) 

< 

> 

SE(A,  B) 

> 

< 

SW{A,B ) 

< 

< 

Table  3:  Direction  Predicates  for  point  objects  in  terms  of  coordinates 


represent  the  relationships  between  two  point  objects  on  X  and  Y  dimensions  respectively.  Ay ,  By  are  the 
y-components  of  A  and  B,  and  AX,BX  are  x-components  of  A  and  B.  Figure  3(a)  illustrates  the  intuition 
of  the  definition.  The  predicates  are  defined  using  direction  equivalence  classes  [22]  by  partitioning  the 
space.  SP(A,B)  means  A  and  B  are  on  the  same  position,  North,  South,  East,  and  West  represent 
exact  directions,  while  NE,  NW,  SE,  and  SW  can  point  to  any  direction  in  their  respective  quadrants. 
This  predicate  set  is  complete.  Figure  3(b)  shows  examples  of  predicates  described  the  directional 
relationships  among  points  A,  B  and  C,  i.e.,  B  is  east  of  A,  A  is  northwest  of  C,  and  C  is  southwest  of 
B. 

The  spatial  constraints  represented  in  terms  of  conjunctions  of  predicates  can  be  maintained  in  two 
graphs:  X-graph  and  Y-graph.  The  nodes  in  both  graphs  represent  the  objects  forming  the  constraints. 
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A 


B 


E 


East(B,  A),  NW(A,  C),  SW(C,  B) 
West(A,  B),  NE(B,  C),  SE(C,  A) 


(a)  Definitions  of  (b)  Examples  of  predicates 

the  predicates 


Figure  3:  Illustration  of  the  predicates 


The  direction  constraints  are  represented  as  directed  edges  in  each  graph  according  to  the  symbol  in  the 
column  2  and  3  of  Table  3.  The  edge  goes  from  the  node  with  smaller  values  to  the  node  with  larger  values. 
If  the  symbol  is  the  two  nodes  are  merged  to  one  node.  Figure  4  shows  the  graph  representation 

c 

V 

D 

X-graph 

Figure  4:  Dimension  graph  for  North(A,B)  A  NW(B,C)  A  SE(B,D) 

for  the  constraint  North(A,B)  A  NW(B,C)  A  SE(B,D).  The  dimension  graph  is  a  union  of  X-graph  and 
Y-graph,  and  is  constructed  according  to  the  definition  of  each  predicate  in  Table  3.  Algorithm  1  is 
the  pseudo-code  of  the  graph  constructing  procedure.  The  input  of  the  algorithm  is  the  conjunctive 
constraint  and  the  output  is  the  dimension  graph.  For  each  predicate  in  the  conjunctionConstraint, 
the  algorithm  invokes  the  sub-function  add_a_predicate_point  to  add  the  predicate  to  the  dimension 
graph.  This  function  calls  function  addNode  to  add  nodes  that  do  not  in  the  dimension  graph  into  the 
graph,  and  calls  function  addEdge  to  add  the  spatial  relationships  between  nodes  into  dimension  graph. 

We  can  easily  summarize  the  computation  complexity  for  this  algorithm.  Let  N  be  the  number  of 
spatial  objects  involved  and  E  be  the  number  of  spatial  predicates  in  the  conjunction. 

•  Time  complexity  =  0(N+E); 

Any  of  the  sub-functions  addNode  and  addEdge  takes  constant  time(0(l)).  findXconstraint 
and  findYconstraint  are  essentially  table  lookup  functions,  which  could  also  be  accomplished  in 
constant  time.  The  whole  algorithm  therefore  have  the  time  bound  of  0(E).  The  generated  graph 
has  at  most  N  nodes  and  E  edges  each  in  X-graph  and  Y-graph. 

•  Space  complexity  0(N+E). 

We  can  process  on  each  dimension  graph  at  a  time,  the  space  requirement  is  also  linear  to  the 
graph  elements. 
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Algorithm  1  Constructing  Dimension  Graph  from  conjunctive  spatial  constraints  for  points:  con¬ 
struct  GraphPoint 

Input:  conjunctionConstraint  is  a  set  of  conjunctive  predicates 
Output:  constraintGraph  consists  of  the  corresponding  X/Y-graphs. 

Graph  constructGraphPoint(Set  of  Predicates  conjunctionConstraint)  { 

Graph  constraintGraph  =  0; 

for  each  entry  p  €  conjunctionConstraint 

add^a.predicate.point^,  ^constraintGraph); 
return  constraintGraph; 

} 

add_a_predicate_point(Predicate  aPredicate,  Graph*  aGraph  )  { 
firstObject  =  getFirstObject(p); 
secondObject  =  getSecondObject(p); 
addNode(firstObject,  secondObject,  aGraph. graphX); 
symbol  =findXconstraint(Table  3,  p); 

addEdge(symbol,  firstObject,  secondObject,  aGraph. graphX); 
addNode(firstObject,  secondObject,  aGraph. graphY); 
symbol  =findYconstraint(Table  3,  p); 

addEdge(symbol,  firstObject,  secondObject,  aGraph. graphY); 

} 

addNode(Node  nl,  Node  n2,  Graph  aGraph)  { 
if  nl  ^  aGraph  Add  nl  to  aGraph; 

if  n2  ^  aGraph  add  n2  to  aGraph; 

} 

addEdge(char  symbol,  Node  nl,  Node  n2,  Graph,  aGraph)  { 
if  (symbol  ==  ’  =  ') 

merge  nodes  nl,  n2  to  one; 
else  if  (symbol  ==  '<’) 

add  directed  edge  of  (nl,  n2)  to  aGraph; 

else 

add  directed  edge  of  (nl,  n2)  to  aGraph; 


2.2  Dimension  Graph  for  Conjunctive  Spatial  Constraints  Among  Intervals 

The  relationships  between  intervals  proposed  by  Allen  [3]  can  be  defined  in  terms  of  the  endpoints  of 
intervals.  Table  4  shows  the  definition  of  the  13  relationships,  where  Ai,A2  and  Bi,B2  represent  the 
start  and  end  points  of  the  intervals  A  and  B  respectively.  Allen’s  symbols  are  used  here. 

The  relationships  represented  in  terms  of  conjunctions  of  predicates  can  be  maintained  in  directed 
graphs,  where  the  nodes  represent  start  or  end  points  of  the  individual  intervals,  and  the  directed  edges 
represent  the  constraints  between  the  two  points.  Each  edge  is  added  to  the  graph  according  to  the 
definition  in  Table  4.  The  edge  point  to  the  nodes  with  larger  values  from  the  nodes  with  smaller 
values.  The  nodes  with  same  values  are  merged  into  one  node.  It  is  worth  noting  that  there  is  an 
intrinsic  constraint  between  the  start  point  and  the  end  point  of  an  individual  interval,  i.e.  start-point 
<  end-point.  Figure  5  shows  the  graph  representation  for  the  constraint  before(S,R )  A  meets(S,  L). 


Predicate  name 

predicates 

point  relationships 

before 

<  (A,B) 

A2  <  B\ 

equal 

=  {A,B ) 

( Ai  =  Bi)  A  (A2  =  B2) 

overlaps 

o(A,B) 

(Ai  <  B\)  A  (^2  >  B\)  A  ( A2  <  B2) 

meets 

m(A ,  B) 

(A2  =  B 1) 

during 

d(A,B) 

( A\  >  Bi)  A  (A2  <  B2) 

starts 

s(A,B) 

( A\  =  B{)  A  (A2  <  B2) 

finishes 

f(A,B) 

( Ai  >  B 1)  A  ( A2  ~  B2) 

after 

bi(A,B ) 

B2  <  A\ 

overlapby 

oi(A,  B) 

( B\  <  A\)  A  ( B2  >  A\)  A  (J?2  <  A2) 

met  by 

mi(A,  B) 

(B2  =  Ai) 

duringby 

di(A,  B) 

(Bi  >  Ai)  A  (B2  <  A2) 

startby 

si(A,  B) 

{B\  =  Ai)  A  (B2  <  A2) 

finishedby 

fi(A,  B) 

(B 1  >  A\)  A  (A2  =  B2) 

Table  4:  Spatial  relationships  for  intervals,  where  <  and  bi  describe  directional  relationships  and  others 
are  topological  relationships 


spatial  constraint 
intrinsic  constraint 


L2 

Before(S,  R)  and  Meets(S,  L) 


SI  S  S2  R1  R  R2 

• - •  • - • 


LI 


L 


L2 


Interval  relationship 


Figure  5:  before(S,  R )  A  meets(S,  L ) 


The  graph  is  constructed  according  to  the  definition  of  each  predicate  in  Table  3.  The  dashed  arrow 
represents  the  intrinsic  constraint  of  the  start  point  and  end  points.  The  pseudo-code  of  the  graph 
constructing  procedure  is  given  as  in  Algorithm  2. 

The  input  of  the  algorithm  is  the  conjunctive  constraint  and  the  output  is  the  corresponding  di¬ 
mension  graph.  For  each  predicate  in  the  conjunctionConstraint,  the  algorithm  calls  the  subfunction 
add  _a_predicate_interval  to  add  the  predicate  to  the  dimension  graph.  This  function  adds  nodes 
that  do  not  in  the  dimension  graph  into  graph,  and  adds  the  spatial  relationships  between  nodes  into 
dimension  graph. 

Let  N  be  the  number  of  objects  (intervals)  and  E  be  the  number  of  interval  predicates  in  Conjunction- 
Constraint.  We  can  summarize  the  complexity  as  follows: 

•  Time  complexity  =  0(E); 

According  to  table  4,  there  are  at  most  three  point  predicates  should  be  satisfied  for  each  interval 
predicate,  and  hence  at  most  three  edges  added  for  each  interval  predicate.  The  time  complexity 
for  this  algorithm  is  therefore  0(E). 

•  Space  complexity  0(N+E). 

The  resulting  constraintGraph  consists  of  at  most  2N  nodes  and  3E  edges.  The  space  requirement 
is  linear  to  the  number  of  nodes  and  edges,  roughly  2N+3E,  which  is  essentially  0(N+E). 
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Algorithm  2  Constructing  Graph  from  conjunctive  constraints  for  intervals:  construct Graphlnterval 

Input:  conjunctionConstraint  is  a  set  of  conjunctive  predicates 
Output:  constraintGraph  is  the  corresponding  graph 

Graph  constructGraphInterval(Set  of  Predicates  conjunctionConstraint)  { 
constraintGraph  =  0; 
for  each  entry  p  €  conjunctionConstraint 

add_a_predicate_interval(p,  ^constraintGraph); 
return  constraintGraph; 


add_a_predicate_interval(Predicate  p  Graph*  aGraph  )  { 
firstObject  =  getFirstObject(p); 
secondObject  =  getSecondObject(p); 

addlntervalNode(firstObject,  secondObject,  constraintGraph); 
pointPredicates  =  convertToPoint(p);  //according  to  Table  4 
for  each  r(k,l)  in  pointPredicates  { 

if  ( k  <  l )  add  directed  edge  (k,l)  to  aGraph; 
if  (k  >  l)  add  directed  edge  ( l,k )  to  aGraph; 
if  ( k  =  l)  merge  node  k  and  l  in  aGraph; 

} 


3  Consistency  checking  for  conjunctive  constraints 

In  the  previous  section,  we  construct  dimension  graph  for  conjunctive  Euclidean  spatial  constraints 
among  point  objects  or  intervals.  In  this  section,  we  will  describe  the  consistency  checking  algorithm 
based  on  the  dimension  graph  representations. 

Let’s  revisit  the  example  given  in  Figure  4.  If  we  add  a  new  constraint  NE(A,  C ),  the  new  constraint 
graph  is  given  in  Figure  6.  A  new  edge  represented  by  dot  line  is  added  into  the  graph.  As  can  be 


A 

B 


X-graph 


B 


Y-graph 


Figure  6:  North(A,B)  A  NW(B,C)  A  SE(B,D)A  NE(A,C) 


seen,  a  cycle  is  constructed  in  the  X  dimension  graph  in  Figure  6.  In  other  words,  the  constraint  of 
North(A,  B)  A  NW ( B ,  C)  A  SE(B,  D)  A  NE(A,  C)  requires  that  the  x-value  of  A  is  smaller  than  the  x- 
value  of  C,  and  at  the  same  time,  the  x-value  of  C  is  smaller  than  the  x-value  of  A.  This  is  a  contradiction, 
which  means  the  constraint  is  inconsistent.  In  general,  we  can  characterize  this  feature  as  Theorem  1. 

Theorem  1  A  conjunctive  constraint  is  consistent  if  and  only  if  there  exists  no  cycle  in  its  corresponding 
dimension  graphs,  i.e.,  the  graphs  are  all  directed  acyclic  graphs(DAG). 
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Proof: 


==>:  Suppose  that  the  dimension  graph  contains  no  cycle,  i.e. ,  the  graph  is  a  directed  acyclic 
graph(dag),  we  can  construct  a  topological  sort  of  the  dag  [7]  using  depth-first  search  algorithm.  A 
topological  sort  of  a  dag  is  a  linear  ordering  of  all  its  nodes  such  that  if  the  dag  contains  an  edge  ( u,v ), 
then  u  appears  before  v  in  the  ordering.  In  other  words,  there  exists  at  least  one  spatial  configuration 
satisfying  the  constraints,  the  constraint  is  consistent. 

<S=:  Suppose  that  the  dimension  graph  contains  a  cycle.  For  the  nodes  involved  in  the  cycle,  there 
exist  no  topological  sort,  i.e.,  no  linear  ordering  of  the  nodes  is  possible.  In  other  words,  we  can  not  give 
a  global  labeling  for  the  nodes  to  meet  the  requirement  that  all  directed  edges  go  from  left  to  right.  No 
spatil  configuration  can  satisfy  the  constrain.  Therefore,  the  constraint  is  inconsistent.  ■ 

3.1  Basic  Algorithm 

We  now  describe  the  algorithm  for  consistency  checking  for  conjunctive  spatial  constraints  among  a  set 
of  points  or  a  set  of  intervals  based  on  dimension  graph  representation.  The  consistency  checking  for 
conjunctive  constraints  contains  two  steps:  1.  Construct  the  corresponding  dimension  graph;  2.  Perform 
cycle  detection  on  each  graph.  The  constraint  is  consistent  if  none  of  the  graph  contains  cycle.  The 
pseudo-code  is  described  as  in  Algorithm  3. 


Algorithm  3  Consistency  checking:  conjunctionConsistencyCheck 

Input:  conjunctionConstraint  is  a  set  of  conjunctive  predicates 
Output:  TRUE  if  consistent,  FALSE  otherwise 

conjunctionConsistencyCheck(Set  of  Predicates  conjunctionConstraint)  { 

Graph  constraintGraph  =  0; 
if  the  constraint  is  among  points 

constraintGraph  =  construct GraphPoint  (conjunctionConstraint); 
else  //the  constraint  is  among  intervals 

constraintGraph  =  construct Graphlnterval  (conjunctionConstraint); 
if  !detectCycle(constraintGraph) 
return  TRUE; 
return  FALSE; 


The  algorithm  first  constructs  the  dimension  graph  by  invoking  subfunction  constructGraphPoint 
or  construct  Graphlnterval  according  to  the  types  of  the  objects.  The  function  detectCycle  performs 
cycle  detection  on  the  corresponding  dimension  graph.  The  algorithm  return  TRUE  if  no  cycle  is 
detected. 

A  nice  property  of  this  algorithm  is  its  efficiency.  The  consistency  checking  is  just  a  graph  cycle 
detection  which  can  easily  be  done  in  linear  time.  As  in  previous  section,  let  N  be  the  number  of  spatial 
objects  and  E  be  the  number  of  spatial  predicates  involved  in  the  conjunctive  constraint. 

•  Time  complexity  =  0(N+E); 

0(N+E)  is  the  time  for  cycle  detection  in  a  directed  graph  with  N  nodes  and  E  edges[7].  As 
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we  explained  in  section  2,  each  of  X-graph  and  Y-graph  for  a  set  of  points  has  at  most  N  nodes 
and  E  edges,  and  the  dimension  graph  for  intervals  has  at  most  2N  nodes  and  3E  edges.  There¬ 
fore,  0(2N+3E)  is  the  upper  bound  time  complexity  for  consistency  checking  for  a  conjunctive 
constraint,  which  is  same  as  0(N+E). 

•  Space  complexity  =(N+E). 

We  can  process  on  each  dimension  graph  at  a  time,  the  space  requirement  is  also  linear  to  the 
graph  elements. 

3.2  Examples 

Consider  the  example  of  spatial  constraints  among  intervals: 

(S  meets  L)  and  (S  is  metby  R). 

The  consistency  checking  for  this  constraint  is  then  accomplished  by  two  steps.  The  function  construct- 
Graphlnterval  is  first  invoked  to  construct  the  corresponding  dimension  graph,  which  is  illustrated  in 
Figure  7  (a).  The  detectCycle  function  checks  the  cycle  in  the  graph.  Since  there  is  no  cycle  found, 

SI  f 

sit 

R2 


(a)  Dimension  graph  (b)  Interval  relationship 


R1  R 


SI 


•  S2 


Figure  7:  Dimension  Graph  for  meets(S,L )  A  metby  (S,R) 


the  constraint  is  consistent.  We  can  construct  an  interval  configuration  which  satisfies  the  constraint  as 
in  Figure  7  (b). 

Assume  a  new  fact  overlaps{L,  R)  is  added  into  the  system.  We  call  the  function  of 
add_a_predicate_interval  to  add  this  new  constraint  to  the  dimension  graph,  resulting  in  a  new  di¬ 
mension  graph  as  in  Figure  8. 


Inconsistent  constraint 


Figure  8:  Dimension  Graph  after  adding  overlaps(L,R) 


The  cycle  detection  is  performed  on  the  new  dimension  graph  to  check  consistency  of  the  status. 
Since  cycles  of  (Rl,  SI,  LI,  Rl)  and  (R2,  LI,  L2,  R2)  are  detected  in  the  figure  as  shown.  The  new 
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constraint  added  is  not  consistent  with  the  system,  therefore,  the  new  constraint  is  inconsistent.  There 
exist  no  interval  configuration  satisfying  this  constraint. 

Since  the  consistency  checking  is  performed  based  on  the  cycle  detection  of  the  corresponding  dimen¬ 
sion  graph,  the  algorithm  can  always  detect  inconsistent  constraints.  As  we  stated  in  Theorem  1,  for 
any  dimension  graph  without  cycle,  there  is  a  consistent  configuration  among  objects.  The  algorithm 
can  guarantee  global  consistency. 

4  Consistency  Checking  for  Conjunctive  Constraints  Among 
MBRs 

In  the  previous  section,  we  discuss  the  dimension  graphs  for  conjunctive  constraints  among  2D  points 
and  ID  intervals,  and  also  describe  the  dimension  graph  based  consistency  checking  algorithm.  In  this 
section,  we  will  extend  our  dimension  graph  based  approach  to  2D  spatial  objects  approximated  by 
MBRs.  Examples  of  2D  spatial  objects  includes  polygon  regions. 

It  is  common  in  spatial  databases  to  approximate  2-D  regions  by  minimum  bounding  rectan- 
gles(MBRs)  which  are  orthogonal  with  respect  to  the  global  coordinate  system.  Figure  9(a)  shows 
a  small  portion  of  the  campus  maps  of  University  of  Minnesota.  Figure  9(b)  replaces  all  the  buildings 
in  Figure  9(a)  by  their  corresponding  MBRs. 


!  Nicholson 
!  Hall 

Pillsbury  ! 

Hall 

V^estbrook!  ; 

Hall  : 

Northrop 

Memorial 

Auditorium 

Johnston  . 

!  Walter  j 
i  Library  | 

!  Tate  Lab  1 
;  of  ! 

;  Physics 

(a)  Campus  map  of  University  of  Minnesota 


(b)  Approximate  buildings  by 
MBRs 


Figure  9:  2D  region  objects  vs.  MBRs 

By  using  MBR  approximation,  we  can  use  two  representative  points,  namely  lower-left  and  upper- 
right  corners  to  determine  the  corresponding  object.  In  the  rest  of  the  paper,  we  use  the  notation  of  An 
and  Aur  to  represent  the  lower-left  and  upper-right  corners  of  the  MBR  for  any  object  A.  The  notations 
of  An  ,,. ,  An  y,  Aur  x  and  Aur  y  are  used  to  represent  the  x  and  y  coordinates  for  lower-left  and  upper-right 
corners  of  MBR  A. 

The  direction  relationship  between  MBRs  can  be  determined  by  the  relationships  between  the  rep¬ 
resentative  points.  Table  5  shows  the  definitions  of  the  direction  predicates  based  on  the  representative 
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points  of  MBR  of  the  objects.  The  first  column  of  the  table  enumerates  the  direction  predicates.  The 


Direction  predicates 

conditions 

SP(A,  B) 

(Aii'X  =  Bu.x)  A  (An  v  =  BumV)  A  (Aur,x  =  Bur,x )  A  ( AUr.y  =  BUr.y ) 

North(A,  B) 

(A u.y  ^  Bur.y )  A  ( Aii  x  ^  Bu  x)  A  (Aur,x  <  Bur,x) 

South(A,  B) 

( AUr.y  5:  ^ll.y)  A  (An.x  ^  Bn.x)  A  ( AUr.x  <  Bur,x ) 

Easi{A ,  B) 

(All'X  ^  Bur,x)  A  ( All  y  ^  Bll.y )  A  ( Aur,y  >  Bur.y) 

W est(A,  B) 

{AUr.x  ^  Bn.x)  A  ( An.y  ^  Bn.y )  A  ( AUr.y  ^  BUr.y) 

NE(A,B) 

( An.x  >  Bur.x)  A  ( An  y  ^  BUr.y ) 

SE(A,  B) 

(An.x  >  BUr.x)  A  ( AUr.y  <  Bn.y ) 

NW(A,  B) 

( AUr.x  ^  Bn.x)  A  ( An  y  >  Bur.y) 

SW(A,B) 

( AUr.x  ^  Bn.x)  A  ( AUr.y  Bn.y) 

Table  5:  Direction  Predicates  for  MBR 


second  column  gives  the  constraints  that  should  be  satisfied  by  their  representative  points. 

For  example,  in  Figure  9(b),  MBR  of  “  Morrill  Hall”  is  north  of  MBR  of  “Tate  Lab”,  and  the  MBR 
of  “Pillsbury  hall”  is  northeast  of  the  MBR  of  “Northrop  Auditorium”.  The  same  relationships  hold 
for  buildings  in  Figure  9(a).  These  are  described  by  predicates  North(Morrill  Hall,  Tate  Lab),  and 
NE(Pillsbury  hall,  Northrop  Auditorium). 

4.1  Dimension  Graphs  for  Conjunctive  Constraints  Among  MBRs 

The  dimension  graph  for  each  direction  predicate  contains  X-graph  and  Y-graph  which  maintain  the 
constraints  that  must  be  satisfied  on  X  and  Y  dimension  respectively.  The  nodes  in  these  graphs  represent 
the  objects  forming  the  constraints.  The  constraints  are  represented  as  directed  edges  according  to  the 
conditions  specified  in  column  2  of  Table  5.  Similar  to  the  X-graph  and  Y-graph  for  points,  an  directed 
edge  goes  from  the  node  with  smaller  values  to  the  node  with  larger  values.  If  the  value  associated  with 
two  nodes  are  same,  the  two  nodes  are  merged  to  one  node.  We  introduce  a  new  ‘thick  arrow’  edge 
here  to  represent  new  relationships  of  <  or  >.  If  the  relationship  between  two  nodes  p  and  q  is  p  <  q, 
a  ‘thick  arrow’  directed  edge  is  added  from  p  to  q.  If  the  relationship  between  p  and  q  is  p  >  q,  a  ‘thick 
arrow’  directed  edge  is  added  from  q  to  p.  The  dashed  arrow  is  used  for  intrinsic  constraint  between  the 
lower-left  point  and  the  upper-right  point  of  a  MBR. 

The  dimension  graph  for  a  conjunctive  constraint  is  constructed  by  adding  the  constraints  specified 
in  each  predicate  in  X-graph  and  Y-graph.  Figure  10(a)  and  (b)  shows  an  example  of  dimension  graph 
for  a  conjunctive  constraint  of  North(A,  B )  A  NE(B,  C )  A  SW{C,  A).  Figure  10(c)  is  a  possible  spatial 
configuration  among  MBR  A,  B  and  C. 

We  now  describe  the  dimension  graph  construction  for  conjunctive  constraints  for  MBRs.  The 
pseudo-code  of  the  graph  constructing  procedure  is  given  as  in  Algorithm  4.  For  each  predicate  in  the 
conjunctionConstraint,  the  algorithm  invokes  the  sub-function  add_a_predicate_MBR  to  add  it  to  the 
dimension  graph.  The  basic  strategy  of  this  function  is  first  to  transform  the  MBR  predicate  to  point 
predicate  set,  and  add  each  point  predicate  into  the  graph.  To  accomplish  this,  addNode  and  addEdge 
are  invoked. 
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(a)  X-graph 


(b)  Y-graph 


(c)  MBR  configuration 


Figure  10:  Example  MBR  configuration  and  its  corresponding  dimension  graphs,  thick  arrow  represents 
<,  thin  arrow  represents  < 

Let  N  be  the  number  of  objects(i.e.  MBRs)  and  E  be  the  number  of  predicates  in  conjunctionCon- 
straint.  We  can  summarize  the  complexity  as  follows: 

•  Time  complexity  =  0(E); 

According  to  table  5,  there  are  at  most  four  point  predicates  should  be  satisfied  for  each  MBR 
direction  predicate,  and  hence  at  most  four  edges  added  for  each  predicate.  The  time  complexity 
for  this  algorithm  is  therefore  0(E). 

•  Space  complexity  =  0(N+E). 

The  resulting  constraintGraph  consists  of  at  most  2N  nodes  and  4E  edges.  The  space  requirement 
is  linear  to  the  number  of  nodes  and  edges,  roughly  2N+4E,  which  is  0(N+E). 

4.2  Consistency  checking  for  conjunctive  constraints  Among  MBRs 

After  constructing  dimension  graphs  for  conjunctive  constraints  among  MBRs,  we  can  check  the  consis¬ 
tency  by  detecting  cycle  in  the  dimension  graphs.  There  are  three  possible  situations  and  corresponding 
results: 

Case  1:  There  exist  no  cycle  (constraint  predicates  are  consistent) 

Case  2:  Every  cycle  consists  of  only  thick  edge  (constraint  predicates  are  consistent) 

Case  3:  At  least  one  cycle  consists  of  non-thick  edges  (constraint  predicates  are  inconsistent) 

Case  1  is  obvious.  No  cycle  means  there  exist  a  consistent  spatial  configuration  among  all  objects. 
Figure  10  is  an  example  of  consistent  constraints  whose  dimension  graphs  contain  no  cycle. 

In  case  2,  since  every  cycle  detected  contains  only  thick  edges.  Recall  that  any  thick  edge  p  — )•  q 
represents  that  the  relationships  between  p  and  q  is  '<’,  i.e. ,  either  p  =  q  or  p  <  q  could  hold.  If  we 
label  all  thick  edges  as  ’=’,  all  the  nodes  involved  in  the  cycle  will  have  same  value.  We  can  then  merge 
this  cycle  with  a  big  node  consisting  of  all  nodes  involved  in  the  cycle.  The  transformed  graph  contains 
no  cycle,  and  the  corresponding  constraint  is  consistent.  Figure  11  shows  such  an  example.  Figure  11 
(a)  is  the  dimension  graph  for  constraint  North(A,  B )  A  South(B,  A),  it  contains  a  cycle  with  only  thick 
edges  in  X-graph.  Figure  11  (b)  is  the  transformed  X-graph  after  merging  each  cycle  to  one  node.  There 
is  no  cycle  in  this  resulting  graph.  This  constraint  is  obviously  consistent.  Figure  11  (c)  is  a  sample 
configuration. 
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Algorithm  4  Constructing  Graph  for  MBR  constraints:  constructGraphMbr 

Input:  ConjunctionConstraint  is  a  set  of  conjunctive  predicates  between  MBRs 

Output:  constraintGraph  consists  of  the  corresponding  X/Y-graphs. 

Graph  constructGraphMbr(Set  of  Predicates  ConjunctionConstraint)  { 

Graph  constraintGraph  =  0; 

for  each  entry  p  €  ConjunctionConstraint 

add_a_predicate  _MBR(p,  ^constraintGraph); 
return  constraintGraph; 

} 

add_a_predicate_MBR(Predicate  aPredicate,  Graph*  aGraph  )  { 
point  Predicates=  findfromTable(a  Predicate); 
for  each  predicate  r(k,  l)  6  pointPredicates 
addNode(/£,  l,  aGraph. graphX); 
addEdge(r(&,  /),  aGraph. graphX); 
addl\lode(&,  l,  aGraph. graphY); 
addEdge(r(&,  /),  aGraph. graphY); 

} 

addEdge(Point  Predicate  pp,  Graph*  aGraph  )  { 
r(k,l)=  getRelationSymbol(pp); 
if  (r  =='<')  add  directed  edge  ( k,l )  to  aGraph; 
if  (r  =='>')  add  directed  edge  ( l,k )  to  aGraph; 
if  (r  =='=')  merge  node  k  and  l  in  aGraph; 
if  (r  =  ='<')  add  thick  directed  edge  (k,  l)  to  aGraph; 
if  (r  =='>')  add  thick  directed  edge  (l,k)  to  aGraph; 

} 


Case  3  is  easy  to  understand.  If  there  is  a  cycle  containing  at  least  one  non-thick  edges,  no  matter 
how  we  label  the  thick  edges  in  the  cycle,  we  can  not  remove  the  cycle.  There  must  be  conflicts  among 
coordinates  of  the  objects  involved  in  the  cycle,  and  hence,  the  constraint  is  inconsistent.  Figure  12 
shows  an  example  of  the  constraints  among  A,  B,  and  C.  The  constraint  is:  North{A,  B )  A  NE(B,  C )  A 
SW(A,  C).  As  we  can  noticed,  there  exist  cycles  in  both  X  and  Y  dimension  graphs  of  the  constraint. 
No  spatial  configuration  of  A,  B,  C  can  satisfy  this  constraint.  The  constraint  is  inconsistent. 

Based  on  the  above  arguments,  we  can  now  describe  the  consistency  checking  algorithm  for  MBRs 
as  in  Algorithm  5.  The  algorithm  first  calls  function  constructGraphMbr  to  construct  the  dimension 


X-graph 


Y -graph 


(a)  Dimension  graph 


•  Au 


A 


B 


(b)  X-graph  after  merging  the  cycle 


(c)  MBR  configuration 


Figure  11:  Dimension  graphs  for  North(A,B)  A  South(B,A ) 


16 


(a)  X-graph 


(b)  Y-graph 


Figure  12:  Dimension  graphs  for  North(A,  B )  A  NE(B,  C )  A  SW{  A,  C ) 


Algorithm  5  Consistency  checking  for  conjunctive  MBR  predicates:  MBRconjunctionConsisten- 
cyCheck 

Input:  conjunctionConstraint  is  a  set  of  conjunctive  predicates 
Output:  TRUE  if  consistent,  FALSE  otherwise 

MBRconjunctionConsistencyCheck(Set  of  Predicates  conjunctionConstraint)  { 

Graph  constraintGraph  =  0; 

constraintGraph  =  construct GraphMbr  (conjunctionConstraint); 
if  IdetectCycle(constraintGraph) 
return  TRUE; 

else  if  the  cycle  contains  thin  edge 
return  FALSE; 
else  return  TRUE; 


graph,  and  then  use  detectCycle  to  detect  cycles  in  the  dimension  graph.  It  returns  TRUE  if  no  cycle 
is  detected.  If  the  detected  cycle  contains  thin  edges,  the  algorithm  returns  FALSE,  otherwise,  returns 
TRUE. 

Let  N  be  the  number  of  objects(i.e.  MBRs)  and  E  be  the  number  of  predicates  in  conjunctionCon¬ 
straint. 

•  Time  complexity  =  0(N+E); 

The  time  complexity  for  cycle  detection  is  0(N+E).  Checking  the  edge  type  in  a  cycle  can  be  done 
easily  by  turning  on  a  flag  if  the  traverse  passes  a  thin  edge.  Therefore,  this  checking  does  not 
require  extra  time  complexity,  the  time  complexity  is  then  0(N+E). 

•  Space  complexity  =  0(N+E). 

The  dimension  graph  consists  of  at  most  2N  nodes  and  4E  edges.  In  order  to  record  the  edge  type, 
each  edge  may  need  an  extra  flag,  which  will  add  another  E  space.  The  total  space  required  is 
roughly  2N+4E+E,  which  is  0(N+E). 

5  Consistency  Checking  for  Constraints  in  General  Format 

In  the  previous  sections,  we  describe  the  dimension  graph  representation  and  the  consistency  checking 
algorithms  for  conjunctive  constraints  among  points,  intervals  and  MBRs.  In  this  section,  we  will  extend 
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the  dimension  graph  based  consistency  checking  algorithm  to  the  constraints  in  general  format. 

5.1  Dimension  Graphs  for  Constraints  in  General  Format 

The  constraints  in  general  format  can  be  transformed  into  Disjunctive  Normal  Form(DNF),  which  is  a 
disjunction  of  conjunctions  where  no  conjunction  contains  a  disjunction.  Each  conjunction  in  the  DNF 
constraint  can  be  represented  by  a  dimension  graph  by  applying  the  algorithm  constructGraphPoint 
or  constructGraphlnterval.  The  dimension  graph  of  the  general  format  constraint  therefore  are 
collection  of  all  the  dimension  graphs  constructed  from  all  its  conjunctions.  The  number  of  graph 
sets(X/Y-graphs  or  interval  graph)  is  the  same  as  the  number  of  conjunctions  in  the  DNF. 


A 

B 


X-graph 


X-graph 


A. 


Y -graph 


(a)  North(A,B)  A  NW(B,C)  A 
SE(B,D)  A  NE(A,C) 


(b)  SE(A,B)  A  NW(B,C)  A  SE(B,D)A 
NE(A,C) 


Figure  13:  X/Y-graphs  for  (North(A,  B)  V  SE{A,  B))  and  NW(B,C)  and  SE{B,D) and  NE(A,C) 

Figure  13  is  an  example  of  the  dimension  graph  maintaining  the  constraint  of  North{A,  B)VSE(A,  B ) 
and  NW(B,C)  and  SE(B,D) and  NE(A,C)  for  point  objects  A,  B,  C  and  D.  The  DNF  format  of  the 
constraint  is:  (North(A,B)  A  NW(B,C)  A  SE(B,D)  A  NE(A,C))  V  (SE(A,B)  A  NW(B,C)  A  SE(B,D)  A 
NE(A,C)).  Figure  13(a)  and  (b)  are  X/Y-graphs  for  the  two  conjunctions,  which  are  ( North(A,B )  A 
NW(B,  C)  A  SE(B,  D )  A  NE(A,  C))  and  (SE(A,  B)  A  NW(B,  C)  A  SE(B,  D )  A  NE(A,  C)). 


(a)  before(S,R)  and  meets(S,  L)  and  overlaps(L,R)  (b)  metby(S,R)  and  meets(S,  L)  and  overlaps(L,R) 

Figure  14:  Constraint  graphs  for  (S  meets  L)  and  (S  before  or  metby  R)  and  (L  overlaps  R) 

Figure  14  is  an  example  of  the  dimension  graph  maintaining  the  constraint  of  ”S  meets  L  and  S 
before  or  metby  R  and  L  overlaps  R”  for  intervals  S,  L  and  R.  The  DNF  format  of  the  constraint  is: 
( before(S ,  R )  A  meets(S ,  L )  A  overlaps(L ,  R ))  V  ( met.by(S ,  R )  A  meets(S ,  L)  A  overlaps(L ,  R )).  The  two 
conjunctions  correspond  to  Figure  14(a)  and  (b)  respectively. 
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5.2  Consistency  Checking  for  General  Constraints 

After  constructing  the  corresponding  dimension  graph  for  the  constraints,  the  cycle  detection  function 
is  performed  on  each  subgraph  representing  a  conjunctive  constraint.  If  each  set  of  subgraph  contains 
cycle,  the  constraint  is  inconsistent.  If  a  subset  of  graphs  contains  cycle,  the  constraint  combinations 
corresponding  to  those  graphs  containing  cycles  are  inconsistent.  The  combinations  corresponding  to  the 
graphs  without  cycles  are  consistent.  Algorithm  6  is  the  pseudo-code  for  consistency  checking  for  general 
format  constraints  among  points  or  intervals.  The  consistency  checking  algorithm  contains  three  steps: 
Normalizing  the  constraints  to  standard  DNF  formats;  Construct  dimension  graphs  for  each  conjunction 
in  DNF;  Perform  cycle  detection  on  each  graph. 


Algorithm  6  Consistency  checking:  consistencyCheck 

Input:  constraint  is  the  constraint  need  to  be  checked 
Output:  TRUE  if  consistent,  FALSE  otherwise 

consistencyCheck(Set  of  Predicates  constraint)  { 

DNF  dnfConstraint=  normalize(constraint); 

Set  of  Predicates  consistentConstraint  =  0; 
for  each  conjunction  p  €  dnfConstraint  { 
if  objects  are  MBRs 

if  MBRconjunctionConsistencyCheck(p) 
add  p  to  consistentConstraint; 
else  if  conjunctionConsistencyCheck(p) 
add  p  to  consistentConstraint; 

} 

if  consistentConstraint  ==  0 
return  FALSE; 
return  TRUE; 

} 


The  subfunction  normalie  preprocesses  the  constraint  and  transform  the  general  format  constraint 
into  its  DNF  format.  Secondly,  check  the  consistency  for  each  conjunction  of  the  DNF  representation 
by  calling  the  Algorithm  conjunctionConsistencyCheck  or  MBRconjunctionConsistencyCheck 
according  to  the  type  of  the  objects.  If  the  conjunction  constraint  is  consistent,  add  the  conjunction 
to  consistentConstraint.  The  final  result  of  consistentConstraint  contains  all  the  consistent  constraint 
combination,  each  of  which  represents  a  global  consistent  constraint. 

In  the  example  of  constraint  among  point  A,  B,  C  and  D  shown  in  Figure  13,  the  dimension  graph 
corresponding  to  conjunction  North(A,B)  A  NW(B,C)  A  SE(B,D)  A  NE(A,C)  contains  a  cycle,  and  hence 
this  conjunction  constraint  is  inconsistent.  The  consistentConstraint  will  only  include  the  consistent 
conjunction,  which  is  SE(A,B)  A  NW(B,C)  A  SE(B,D)A  NE( A, C) (Figure  13(b)).  Similarly,  in  the 
example  for  intervals  illustrated  in  Figure  14,  the  conjunction  depicted  in  Figure  14(b)  contains  cycle 
and  hence  the  corresponding  conjunction  is  inconsistent.  Figure  15  shows  the  only  consistent  conjunction 
and  one  of  its  possible  interval  configurations. 
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Figure  15:  S  meets  L  and  S  before  R  and  L  overlaps  R) 


Complexity  Analysis 

Similarly,  let  N  be  the  number  of  spatial  objects  and  E  be  the  number  of  spatial  predicates  being  checked 
for  consistency,  and  COR  be  the  number  of  all  possible  combinations  of  disjunctive  predicates,  i.e. ,  the 
number  of  conjunctions  in  DNF. 

•  Time  complexity  =  0(N+E)*0(C0R); 

The  consistency  checking  algorithm  for  general  format  constraints  first  generates  the  DNF  for  the 
constraints,  and  then  call  the  conjunctionConsistencyCheck  algorithm  for  each  conjunction. 
0(N+E)  is  the  time  for  consistency  checking  for  a  conjunction  constraint.  Therefore,  the  total 
time  bound  is  0(N+E)*0(C0R). 

•  Space  complexity  =  0(N+E). 

COR  does  not  contribute  to  space  factor,  since  one  may  process  one  graph  at  a  time. 

6  Conclusions  and  Future  Work 

In  this  paper,  We  propose  a  new  strategy  to  process  consistency  checking  for  Euclidean  spatial  con¬ 
straints  among  objects.  We  use  a  geometric  approach  by  incorporating  the  spatial  domain  information 
in  consistency  checking.  Dimension  graphs  are  proposed  to  maintain  the  spatial  constraints  among  ob¬ 
jects.  Each  conjunctive  constraint  is  projected  on  both  dimensions,  and  a  dimension  graph  is  constructed 
for  the  constraint  on  each  dimension.  The  spatial  constraints  in  general  format  are  maintained  in  a  set 
of  dimension  graph  constructed  from  its  conjunctions.  By  using  this  framework,  constraint  consistency 
checking  is  then  transformed  to  a  graph  cycle  detection  problem  on  its  dimension  graph.  The  cycle 
detection  could  be  solved  efficiently  with  0(N+E)  time  as  well  as  space  complexity,  where  N  is  the  num¬ 
ber  of  spatial  objects,  and  E  is  the  number  of  spatial  predicates  in  the  constraint.  Recall  that  Allen’s 
algorithm  has  a  time  and  space  complexity  of  0(N2).  The  proposed  approach  to  consistency  checking 
for  spatial  constraints  is  faster  when  the  number  of  predicates  is  much  smaller  than  N 2  and  there  are 
few  disjunctions  in  the  spatial  constraint.  The  dimension  graph  and  consistency  checking  algorithm  can 
be  used  for  points,  and  MBRs  in  2  dimensional  space.  Since  the  algorithm  returns  TRUE  if  and  only  if 
the  dimension  graph  of  at  least  one  conjunction  contains  no  cycle,  which  means  there  exist  at  least  one 
consistent  overall  constraint.  The  algorithm  guarantees  the  global  consistency. 

In  future  work,  we  would  like  to  explore  the  consistency  checking  among  mixed  types  of  object, 
e.g.  consistency  checking  for  constraints  between  point  and  intervals.  We  would  also  like  to  apply  the 
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dimension  graph  based  cycle  detection  algorithm  to  the  application  of  image  similarity  retrieval  based 

on  the  similarity  of  spatial  configuration. 
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